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Abstract 

We show that the Frobenius group T13 = Z\% x Z3 is a suitable family symmetry group, to 
study neutrino oscillations. Our approach is to catalog all possibilities within an effective field 
theory approach, assuming only SU(2)xU(l), supplemented by family symmetry. We will use 
tribimaximal mixing as a guide to place a constraint on the otherwise various possibilities. This 
leads to an exact fit between the neutrino and charged lepton sector. Such a fit has not been 
achieved with any other group so far. The results of this paper may then be useful in future 
studies on the compatibility of this Frobenius group with other models and mechanisms. 
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1 Introduction 

According to the Standard Model, neutrinos do not have a mass. However, experiments have 
shown that neutrinos oscillate, and in order to do so, there must be a difference in their masses, 
so that at least two of them are massive. Evidence that neutrinos undergo oscillations was 
first suggested in the 1960's by R. Davies where a lack of observed solar electron neutrinos 
and atmospheric muon neutrinos compared to the expected amount was measured. Several 
experiments have since then confirmed this deficit, and in 1998, the existence of neutrino masses 
was manifested from observed atmospheric neutrino oscillations [lj. 

The existence of a family symmetry group F (also called horizontal symmetry) was pro- 
posed long ago p|. Under F, the 3 quark and lepton generations transform into each other. 
Considerable effort has been put into deriving the mixing and mass matrices for neutrinos, us- 
ing non-abelian finite groups. The tetrahedral group A4 has been studied extensively [SHE]. So 
far, this group seems to present the best fit to describe neutrino oscillations. However, other 
candidates can be found. The Frobenius groups, as sugroups of SU(3), seem to provide interest- 
ing family symmetries. We will focus on the Frobenius group X13 = Z13 x Z3 in this paper. As 
usual, the charged lepton fields (left and right handed), the neutrino fields, and various Higgs 
fields are to be assigned to various representations of T13. Interestingly, T13 has 4 different 
3-dimensional representations, which thus allow for a host of possibilities for model building. It 
has been shown in [9] that T13 is suitable to produce tribimaximal mixing. Instead of making a 
more-or-less random choice, we catalog and study the various possibilities, using the tribimaxi- 
mal mixing matrix as a guide. We have subsequently learned that the group T13 has also been 
discussed in [ID]. This paper makes a particular choice of representations and also proposes a 
theory of decaying dark matter based on this group. We also note in this context that another 
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Frobenius group, T7 = x Z3, has been discussed in and |12j in connection with neutrino 
mixing. 

We will start by a brief review on neutrino mixing and then discuss the group theory of 
the Frobenius group T13. After describing how to extract the charged lepton and neutrino 
mass matrices from the Lagrangian, we will proceed to the model building, where we assign 
different representations to the respective fields involved, and work out some interesting cases. 
We show, that this Frobenius group, with tribimaximal mixing as a guideline, has an even better 
fit between the charged lepton sector and the neutrino sector, than the tetrahedral group. 



2 Neutrino mixing 

Neutrino oscillations require a mixing matrix that transforms the mass eigenstates into the flavor 
eigenstates as follows: 



This mixing matrix is given by: 
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The three k phases can be absorbed into the other phases, by rephasing the neutrino field |13j . 
For Majorana neutrinos, besides the Dirac phase 5, there are two extra Majorana phases, $1 
and < I>2- These are relative phases among the Majorana masses. The Majorana phases are not 
observable in neutrino oscillations, as is well known. 

Experimental measurements are well described by the tribimaximal mixing matrix [Tl] , which 
we will use as a guide for model building in this paper: 
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3 The family group 

The Frobenius group T13 is a subgroup of SU(3) generated by 2 elements of SU(3), a and b, 
such that a 13 = I, b 3 = I and ba = a 3 b. All elements g of X13 can be written as g = a m b n with 
< m < 12 and < n < 2. Thus, the group has 13 • 3 = 39 elements and 7 conjugacy classes 
given by: 

Ci : {e} 

Cg 5 : {b, ba, ba 2 , ... ,ba n ,ba 12 } 
C 13 : {b ,b a,b a , ... b a ,b a ) 
C 3l : {a, a 3 , a 9 } 
C Sl :{a 4 ,a 10 ,a 12 } 
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C 32 : {a 2 , a 5 , a 6 } 

C3 2 :{a 7 ,a 8 ,a n } (4) 

From Sjd? = 39, where di are the dimensions of the 7 irreducible representations (IR) and from 
Ejmj = 7, where mi is the number of IR with dimension di, we see that the group consists of 
three 1 dimensional IR's 1, 1' and 1' and four 3 dimensional IR's, 3i,3i,32 and 32- Explicitly, 
for 1' and 1', a is represented respectively by uj and oj 2 and b by 1, where u> = e~3~ is the cube 
root of the identity. In the four 3 dimensional IR's, a is represented by 



oi = p 3 (5) 
V p 9 J 

( p 2 o \ 

a 2 = p 6 ] (6) 
and the corresponding complex conjugates: 
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where p = e^3~ is the IS*' 1 root of the identity. 

Another way of defining the group T\% = Z13 x Z3 is to demand that for any element b of Z3 
on any element a of [15], we have: 

bah' 1 = a r (10) 
with r some integer. Letting the element of Z3 act on the element of Z13 three times we get: 

b 3 ab- 3 = cT 3 =a (11) 

So we have r 3 = 1 + 13n, with n some integer. Clearly, r = 3 and n = 2 solves this equation, so 
that: 

bab" 1 = a 3 (12) 

(r 7^ 1 is indeed what the symbol x means). The Frobenius group T13 thereby has the presen- 
tation: 

(a, b\a 13 = b 3 = I, bab~ l = a 3 ) (13) 



To find the Kronecker products and the Clebsch-Gordan decompositions of the IR's, it is 
convenient to use the fact, that the Frobenius group is a subgroup of SU(3), and exploit the 
well-known properties of SU(3). Denote the 3-component vector of the defining representation 
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3 by ip 1 *, with \x = 1,2,3, and the vector of the conjugate representation 3 by ip^ = {ip 11 )* . In 
SU(3), we have 

3^3 = 306 (14) 

as indicated by 

i/j^ u ~ ^ v ~ ^"e^ + ^ } (15) 
Under the action of a generator c, ip^ v transforms as follows: 

^ = c^ c \r p (i6) 



To find the decomposition of 3i ® 3i, we use the generators a\ and b to determine how ip^v} 
decomposes. 

Consider the diagonal elements ip^- 11 ^ , tp^ 22 } , ip^ 33 \ Under b, they transform as 

^{11} ^ ^{22} ^ ^{33} ^ ^{11} (17) 

while under a±, they transform as 

^{n} V {22} ->pV {22} , V {33} ^ pV 33} (18) 

Thus they form a 32 representation according to Q: 

^,{11} 

V^ 22 > I (19) 



^{33} 



3 2 



Next, consider the off-diagonal elements of the symmetric representation ip^ u ^ . Under b they 
transform as 

^{12} ^ ^{23} ^ ^{31} ^ ^{12} (2Q) 

Under a±, they transform as 

^,{12} p 4^{12 }) ^{23} ^ p 12 v ,{23} ! ^{31} ^ ^{31} (21) 
thus forming the 3i representation according to Q: 

^,{23} V 

V^ 31 * (22) 
^ {12} / Sl 
We thus obtain 

3i ® 3i = 3i ® 3i ® 3 2 (23) 

The decompositions of 3 2 <8) 3 2 and 3i ® 3 2 are found in the same way, using ai and a2 appro- 
priately. The result is displayed in Table [TJ 

In SU(3) we have the product between a 3 and a 3 representation: 

3 <g> 3 = 1 ® 8 (24) 

where ip^ip u ~ ipv is decomposed into a trace tpfi that forms the singlet and a traceless ^ that 
forms the adjoint 8 representation. 
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For 3i (g) 3i, the diagonal elements of V>v transform under b as f/'i — > — > — )• We 
have the two possibilities 

= ^+a^l + w 2 V>I (25) 

= Vi + w 2 V>| + wV| (26) 

Under b, these transform as 

Vl + ojtpl + wV} = u 2 4> (27) 

-»• -0| + w 2 V| + w$ = (28) 

Thus, we choose (j) to transform as 1' and <j) to transform as 1'. 

The off-diagonal terms split into two sets. They transform under b in the same obvious way 
as described in ( 20 ) . Under a\ and a\ they transform as 

^^p 11 ^ (29) 

and 

ift^fift, Vl^/A^L ^pVs (30) 



giving a 32 and a 32 representation: 



4 | (3D 

^ 3 J (32) 

^ / 3 2 



The rest of the decompositions can then be found in the same way using the appropriate 
representations of the generators to investigate how the elements transform. The products 
between the ID IR's are found from seeing that under 6, (ft — > uxf) and 4> — > lo 2 (J), so that for 
example cpcj) — > uxpuj 2 ^) = <P4> transforms as a 1. The Kronecker products are listed in table [l] 
and the Clebsch-Gordan decompositions in table ^1 where we have dropped the {..} on tp. 
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Table 1: Z13 » Z3 Kronecker products 
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TaWe 2: Z13 xi Z3 Clebsch-Gordan decompositions. 



4 Obtaining the mixing matrix 

We assume in this paper, that the neutrinos are Majorana type. To acquire mass, the neutrinos 
couple to one or more complex Higgs fields through Yukawa couplings. We have the following 
mass terms for the neutrinos v and charged leptons I in the Lagrangian: 

Cmass = -lLa(Mi) a pl R p - ^v[ a {M v ) a? V L ^ + h.C (33) 

where a, (3 denote the flavour indices e, fi, r. The masses are then obtained by finding the 
eigenvalues of Mi and M v through diagonalization. 



V\M X V R = A = (34) 
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U^M U U V = D v = y m 2 | (35) 

The left handed charged lepton mass basis is given by: 

ha = (U L ) ai l Li (36) 



and the neutrino mass basis by: 



with i = 1, 2, 3 as in Q, so that: 



v a = (U u ) ai Ui (37) 



Thus, the lepton mixing matrix is given by U = Jj\U u . 

In this paper, we adopt an effective field theory approach within a minimalist framework, 
assuming only SU(2) <S> U(l) as in [TJ. The charged lepton masses are generated by the dimension- 
4 operator : 

O a = <\>H c i) (39) 

with <j) the Higgs doublet, ijj the lepton doublet and l c the charge conjugated lepton. Neutrino 
masses are generated by the dimension-5 operator |16j : 

5 = (6t-2^)C(6t 2 V) (40) 

Within the Standard Model, this dimension 5 operator is the lowest dimension operator which 
can give masses to the neutrinos. Here the Higgs doublets £i and £2 may or may not be identical 
and there can be as many as needed to carry out the model building [T7J- We will assume that 
they are distinct from the Higgs doublet </>. When the different (or equal) Higgs doublets £1 
and £2 acquire a vacuum expectation value (vev) for their lower electrically neutral components, 
these will give masses to the neutrinos through a symmetric mass matrix. 

In this paper, we do not write down the Higgs potential explicitly and minimize it. With 
several Higgs doublets, the potential will contain a multitude of parameters, and we will not 
learn much by explicitly showing how the vev's of the various Higgs fields are aligned in various 
regions of the parameter space. Instead, we systematically list the various possibilities and an- 
alyze the consequence for the neutrino mixing matrix. 

We assume in the following, that there is no coupling between the Higgs field in the charged 
lepton sector and the Higgs field(s) in the neutrino sector. In other words, we have not adressed 
the sequestering problem [18] |19j . Furthermore, we assume CP conservation, so that M v is real. 



5 Model building 

Given the review of group theory in section 3, we are now ready to build a model for neutrino 
mixing based on the family symmetry group F = Z13 x Z3. As in the standard procedure, we 
assign the lepton doublet fields ip a (a = 1,2,3), the right handed lepton singlet fields l c a , and 
the Higgs fields to various representations of X13, and write down the various allowed forms of 
the operators O4 and O5. After the Higgs fields are allowed to acquire their vev, the resulting 
charged lepton mass matrix M\ and neutrino mass matrix M v are diagonalized, whereby the 
mixing matrix is obtained as outlined earlier. In order to constrain our systematic search, we 
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will use the following assumptions during the model building: 

i) We will align the vev of the Higgs field in the dimension-4 operator in the (1,1, 1) T direc- 
tion. 

ii) We will assume no coupling between the Higgs fields in the charged lepton sector and the 
neutrino sector. 

iii) Since tribimaximal mixing gives a fairly good description of the experimental mixing matrix, 
at least as a first approximation, we will take it as a phenomenological guidepost. 

We have in mind a program in which we could eventually perturb the possible schemes, we 
obtain in this way, so as to study deviations from tribimaximal mixing |19j |20j . 

In the charged lepton sector, the operator O4 = (f>H c ip gives a priori 5 3 = 125 ways of placing 
the three fields in different representations, but some of these are not allowed by the group 
theory. Also, two cases, where the representations are the conjugate of each other, will give 
the same result, so the number of different forms O4 can acquire, is down to 27, see Table |3| 
For each of these cases, we then consider the various possibilities for O5 assigning £1 and £2 
to various representations of F. We then use tribimaximal mixing as a constraint. It turns 
out that the cases 1, 11, 22 and 25 will lead to the same neutrino mass matrix. Another mass 
matrix can be found from the cases 23 and 24. Furthermore, the rest of the cases turn out to 
create degenerate charged lepton mass matrices, which are therefore not of phenomenological 
importance. Therefore, we can cover all possible results obtained from this group, by considering 
only four of the mentioned cases, where the lepton doublet is in two different representations. 
This leads to only two possible neutrino mass matrices. These cases will now be studied. In the 
following, we will not distinguish between upper and lower indices. 
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TaMe 5: Possible representations for the 3 fields. 



6 Case 1 

The charge conjugated lepton, l c , and the lepton doublet, ip, both in 3i can be combined in two 

ways, 

3i: {^3,^1,^2} 

3i: {^2,^3,^1} (41) 
With ft in 3i, we get the invariant terms in the dimension 4 operator of the Lagrangian: 

hi(4l C 2^3 + + 4?M + h 2 (4>[^2 + 4^3 + 4Wi) (42) 

/ o htfl h 2 ft 2 \ / Vi 

h 2 4 Mi ^ 
\ hx4 h 2 (fti J \ih 

Upon spontaneous symmetry breaking, we let the Higgs field acquire a vev in the (1, 1, 1) T 
direction, and the mass matrix becomes: 

/ hi h 2 \ 
Mi = v I h 2 h ± (43) 
\h ± h 2 / 

To diagonalize M/, we multiply by its hermitian conjugate: 



Cl) ^23 ^3 
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/ I h | 2 + | h 2 ? h\h 2 h x hl 

MiM} = M\M l = hih* 2 | hi | 2 + | h 2 | 2 h\h 2 j (44) 

V h\h 2 h x h\ | hx | 2 + | h 2 | 2 

This is a circulant matrix, that can be diagonalized by the following unitary matrix: 

c/t = _ i , i ^ j m 

We now look at the dimension-5 operator and ask what Higgs field vev's will lead to tribi- 
maximal mixing Q The Higgs fields ft and & in @) can each be any of the following 5 different 
Higgs in the representations: 

X ~ 1, T] ~ 3i, p~32, k~3i, C ~ 32 (46) 

We have a priori 5 + 4 + 3 + 2 + 1 = 15 combinations of these in the dimension-5 operator. 
In Case 1, we have the ip in the 3i representation. From Table [TJ we can already exclude the 
combinations (x^) 2 > i 7 !^) 2 -, (xVO^VO an d (xV'XWO- The possible invariants from combining 
two lepton doublets with two Higgs in the above mentioned representations and contracting the 
indices properly, can now be cataloged. 



a) (W) 2 : 

b) (ra/0 2 : 



c) (CV) 2 : 

d) (xV>)(#) : 

e) (xV)(C# 

f) (^(pi/j): 



g) (#)(/#): 

h) M)(CV): 



fl(P3^2Pli>2 + Pl1p3P2lp3 + p2lplP3^l) (47) 

/ 2 («l^lKl'i/'l + K 2 1p 2 K 2 1p 2 + K 3 lf) 3 K 3 i) 3 ) 
+ f3(nitplK 2 7p 2 + K3lp 3 Ki1pl + K 2 1p 2 K 3 1p 3 ) 

+ f4,(Kllp2K>2lpl + K 2 1p 3 K 3 7p 2 + K^l^l^) (48) 

/s(Ci^iCi^2 + (2^2(2^3 + CsV'sCa^i) (49) 

+ }7x{^im^2 + ^2Vi^3 + ip3m^i) (so) 

/sXO/'lCl^l + ^2^2 + ^3^3 ) (51) 
h{m^3P3'>p3 + V3lplPltpl + m4>2P2lp2) 

/10 (^3^2/53^3 + mVw^i + m^ipzfa) 
+ /n {mfap&^i + vifapiipa + mi>2p2^i) (52) 

/12 (vitpi^ipi + m^^ipi + 773^3^1^3) (53) 

/l3 (^3^2(1^2 + 771^2^3 + ^2^1(3^1) 



In order to find the different invariant terms that enter the Lagrangian through this operator, it is necessary 
to keep in mind, that since the neutrinos are Majorana, we have here used the 2 component spinor notation, so 
that each lepton and scalar carries the Lorentz spinor indices a = 1,2. For example, ^1/1 = t^ a ipa ~ tapip 01 ^ 
|21j . We also have to take into account that the leptons are Grassman numbers, and thereby anticommute. 
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i) {pip)(Kifj): 



j) (WO(C# 



k) (K^)m- 



+ h^{j]i^iC^2 + m^Ci^z + m^3(2^i) (54) 

+ fl7(pi1p3K\1p2 + P2^\K-2^3 + P3^2^\) (55) 
/l8 002^3(3*02 + P3^lCl^3 + Pl^Ql^l) 

+ /ig(P3'03CiV'i + Pi^ibfa + ^2^2(3^3) (56) 

f2Q{^2Czi>2 + K2Tp3(l1p3 + ^3^1(2^1) (57) 

Of course, in building a specific model, we are free to omit one or more of these Higgs fields, 
narrowing down the possibilities further. There should to begin with be 5 different terms in b), 
since 3i © 3i = 1 © 1' © 1' © 32 © 32 multiplied by itself creates invariants from 1 © 1, 1' © 1', 
32 © 32 and two charge conjugates. However, when multiplying out the terms and seeing that 
oj + uj 2 = —1, we get the three listed terms. 

In Case 1, we first consider including only the Higgs fields r\ and k, to which we assign the 
vev's: 

i) 



/ON 








j i («) = «1 j 
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V o J 







(V) = v X I 1 J , (k) = u x y I (58) 

so that we get contributions from invariants in b) and g). We could also add % so we have the 
vev's: 











j , («) = m j 
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(X) = <»7> = 1 , («) = «i (59) 



(with 7) and k having the same vev's as in (J58j) ) . This would give another contribution from d) . 

Upon spontaneous symmetry breaking, the neutral lower components of the SU(2) Higgs 
doublets can acquire these vev's and so will connect with the neutrinos giving masses to these. 
The charge conjugation operator, together with the neutrinos being Grassman, will then allow 
us to combine some of the coupling constants in b), d), f), h), i) and j), contributing to the 
off diagonal elements of the neutrino mass matrix, thus making this symmetric as it should be. 
In the following equations, we will, in some cases, combine f a and f a +i into a new coupling 
constant, which we call f a . The terms contributing to the neutrino mass matrix in i) are then 
given by: 

hu\{l>\V\ + V 3 V 3 ) + f^VxVj, + h2V\U\V 2 V2 (60) 

and for (x) getting a vev in ii), we also get the contribution: 

f' & XVlV\V3 (61) 

This leads to the possible mass matrices: 
i) 

hu\ o r z u\ 

M u = | f 12 v lUl | (62) 

f'A o f 2 u\ 
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and 



f2u\ &u\ + 

M„ I favxui 

f 3 u\ + ftxvi hu\ 

with and without the Higgs \ respectively, both having the form: 



(63) 



which is diagonalized by the unitary matrix: 



We can now find the lepton mixing matrix: 
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(64) 



(65) 



U 



u\u v 




(66) 



This is an interesting result. The neutrino mass matrix in (|64|), which leads to tribimaximal 
mixing in case 1, is exactly achieved, by choosing the right vev's in the neutrino sector. This was 
not possible in [7] for the tetrahedral group. In this paper, when choosing vev's for the Higgs 



fields in the neutrino sector, the exact form of the mass matrix in (64) could not be achieved. 
The mass matrix had to be perturbed in the following way: 



a + e (3 
M„ I 7 

(3 a - e 



(67) 



which did not allow for exact tribimaximal mixing. It was instead assumed, that e being small, 
could be seen as a perturbation, and thereby only lead to a small deviation from tribimaximal. 
If tribimaximal mixing is the best fit for lepton mixing, we already here have evidence, that the 
Frobenius group can be used to describe neutrino mixing. In order to determine which of the 
two groups, the tetrahedral or the Frobenius, is the best candidate to describe neutrino mixing, 
it will be necessary to test the Frobenius group further, using more constraints. 



7 Case 23: 

The lepton doublets are in the 3i representation, the charged lepton singlets in the 3i represen- 
tation, and the Higgs bosons (ft in the 1, 1' and 1' representations. 3i and 3i can be combined 
to form the 3 1 dimensional representations: 

1 ~ + ^2 + V>3^ 

1' ~ fall + Ulfal% + oj 2 ^11 

r ~ mi + uPfaq + u^ii (68) 
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Combined with the <f>\ we get the following contributions: 

hi4>\ + h 2 uJ<t>2 + h 3 uj 2 (j)l 

(i c i,i c 2 ,i c 3) | o mI+mS + mS ° 

+ h 2 uj 2 (f)\ + h 3 uj(f)l ) 

1P2 (69) 

^3 / 

When the Higgs field $ again acquires a vev in the (1, 1, 1) T direction, the charged lepton 
mass matrix is already diagonal. To get tribimaximal lepton mixing, we then need the neu- 
trino mixing matrix to be tribimaximal. We use the invariants as listed for case 1 above for the 
if> in the 3i representation to achieve this by letting the Higgs bosons acquire the following vev's: 

i) 

(V) = vi | 1 j , (p) = v 2 ^ 

(70) 

This will give contributions from b), f), g) and i). Letting x acquire a vev in addition to these 
could also be a possibility: 

ii) 

°\ f 1 

(\) - .»-. (//) - | 1 I, (p) = v 2 I 





(71) 

thereby gaining a contribution from d). Writing up these terms in the neutrino mass matrix, 
we find: 



i) 

and 

ii) 



f 2 u\ + faviv 2 f&l f z u\ 

M v = ( f z u\ j 2 u\ + f 12 v lUl f z u\ + f w v 2Ul \ (72) 

f!i U l f'z U l + fl6 v 2Ul + f\ 2 V\Ux 



hu\ + hv\v 2 fzii\ + f' & xvi f^ul + f & xvi 

M» = | &u\ + &x Vl f 2 u\ + fav lUl f z u\ + f[ 6 v 2Ul I (73) 

f3 u l + f'% xv l f3 u l + fl6 v 2 u l Hu\ + f\ 2 V\U X 

leading to a matrix with p — r symmetry: 

( 1 P P\ 

M u = \ ft a 8 (74) 
V P 5 a J 

which is diagonalized by the tribimaximal mixing matrix, if we impose without justification, 
5 = 7 + f3 — a. 
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8 Case 25: 

We now have the l c in the three ID representations, the lepton doublet in the 32 representation 
and the (ft in the 32 representation. The lepton doublet and the Higgs then join to form the 
three 1 dimensional representations: 

i ~ 4v>i + 4>\i>2 + 4^3 

v ~ 4v>i + w 2 4^2 + ^4^3 (75) 

which results in invariants when connected with 1 ~ I \ , 1' ~ Z3 and 1' ~ Z2 respectively. The 
mass terms of the charged leptons in the Lagrangian are then (using the freedom to incorporate 
phases into the coupling constants hi): 

hilt(<&i/)i + 4V>2 + 4^3) + hsuPlittlifa. + w4V>2 + w 2 4^s) 

+ h2U)l^(<f>\lf>l + W 2 4^2 + U^lh) = ^1^(4^1 + 4^2 + 4^3) 
+ h 2 l C 2 {0j4>\'4)l + 4^2 + W 2 4^3) + /l3^(w 2 4V'l + 4^2 + ^4^3) 

hi \ / 1 1 1 

( 11 l c 2 1% ) [ h 2 
h 3 





4 I I ^2 I (76) 

Upon spontaneous symmetry breaking, again with the Higgs vev in the (1, 1, 1) T direction, the 
mass matrix is diagonalized by the same unitary matrix as in case 1: 



1 1 1 

Proceeding to the neutrino sector, we now have tp in the 32 representation, leading to different 



Ul = ^=\ w 1 ^ ] (77) 

CO 2 1 UJ 



invariants in the dimension-5 operator than for case 1 and 3. The Higgs fields £1 and £2 hi (40) 
can each still be any of the following 5 Higgs with the assigned representations: 

J7~3i, p~3 2 , k~3i, C~ 3 2 (78) 

Using table[l]once again, we can omit the combinations (xV0 2 > (/°^) 2 > (xVOOCVO an d (xV0( K ^)- 
This leaves the possibilities: 

a) (t?V0 2 : 



b) (ra/0 2 : 

c) (CV) 2 : 



d) (xV0W) : 



himi'im^ + m^vi^ + m^mtpi) (79) 

f 2 {Kl^lK 2 lj)i + K2^2^2 + K 3 1p 3 K,l1p3) (80) 

/3(ClV , lClV'l + C2V'2C2V ; 2 + CsV^^) 
+ /4(ClV ; lC2V'2 + C3^3Cl*01 + (2^2(3^3) 

+ MClV^V'l + CsV'lClV'S + C2'03C3V'2) (81) 

kxi^im^i + ^2Vi^2 + ^3*72^3 ) (82) 
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e) (xV)(P# 



f) (th/0(p# 



g) 

h) (^)(C^) = 

i) (pil>)(Kip): 



j) (pYO(C# 
k) (k^)(C# 



+ fsXi^lPS^ + ^2/91^3 + 1p3P2^l) (83) 

h{m^2Pz^2 + m^pi^ + m^iP2^i) 
f\o{m'4 ) 2p2^ + vi^pstpi + m^ipifa) 
+ /n {m^P2^2 + m^ip^z + f]2i"iP\^\) (84) 

712(^3^1^2^3 + vi^Ksipi + m^Ki^) 

+ /l3(%V'3K2'01 + Vl^l^2 + mip2Kltp3) (85) 

/i4(%V'3C2V'3 + m^iC^i + m^Ci^) (86) 

fl5(p3lp2K 2 1p3 + Pl^3/«3^1 + P2^\K\^2) 
+ fw(P2lp2Kllpl + P3^3^2 + Pl^lK 3 ^3) 
fn(P2^3H2lp3 + P3^l^l + Pl1p2Kl1p 2 ) (87) 



hs{P2^2C,3^2 + ^3^3(1^3 + ^1^1(2^1 



/l9(«lV'lC2^3 + ^2^3^! + ^3^3(1^2) 
+ /20(«2^lC3^2 + K 3 ^ 2 Cl^3 + KlV^V'l) (89) 



Evidently, the coefficients here are not the same as for ?/> in the 3i representation. In this case, 
the possible Higgs field vev's needed to eventually get tribimaximal mixing are numerous. We 
list them all in the following: 

i) 

<P> = ^2 ^ 1 j , <C) = «2 ^ j (90) 

leading to contributions from c) and j). 
ii) 

( X ) =x, (p) = v 2 ^ 1 j , (C) = n 2 ^ j (91) 

which gives the same contributions as in i) with one extra from e). 
iii) 

<K> = U! ^ 1 j , (C) = «2 ^ 1 j (92) 

with b), c) and k) contributing to the mass matrix, 
iv) 

( X ) = x, (k) = m ^ 1 j , (C) = us ^ 1 I (93) 
leading to the same result as in iii) with no extra contribution. 
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v) 



(/>) - r-2 I I , {k)= Ui I 1 | (9 1.) 



giving terms to the mass matrix from b) and i). 
vi) 

(V) = vi ^ 1 j , (C) = n 2 ^ | (95) 

contributing with terms from a), c) and h). 

When the neutral component of the Higgs bosons acquire these vev's we can then combine 
the coupling constants as we did for the invariants with ijj ~ 3i. The possible mass matrices 
from these couplings and vevs will then become: 



i) 



ii) 



hi), iv) 



v) 



fou\ fou\ 
M v = | f 18 v 2 u 2 | (96) 

fi4 o hul 

foul f A ul + fljxv 2 

M v = | f ls v 2 u 2 ) (97) 

fou\ /{ 9 uiu 2 

M„ = | fou\ + fou 2 2 J (98) 
/( 9 uiu 2 / 2 u? 

M„ = | fouj + f 17 V2U! I (99) 
f[ 5 v 2 ui / 2 «i 



VI 



/ 3 «1 / lWl 2 + 

M v = | /i4«i«2 I (100) 

fivf + fiul / 3 ul 

all with the same form as in case 1: 



a 





P 





7 





P 





a 



M v = 7 (101) 
\ p a / 

and as in case 1, tribimaximal mixing is possible from this mass matrix. 

9 Case 24: 

This case gives the same result as case 23 in the charged lepton sector. We therefore have the 
following contributions to the charged lepton mass matrix: 



Ol) ^21 ^3) 



( hi<j>\ + h 2 Lu4>l + h 3 uj 2 (j)l \ 

hKf>{ + h 2 (f>\ + M3 

hi(f>\ + /i 2 w 2 4 + h 3 u(f>l J 
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(102) 

where the charged lepton mass matrix is already diagonal upon alignment of the $ vev in the 
(1, 1, 1) T direction, as in case 23. We therefore want a tribimaximal neutrino mixing matrix 
once again, now with the ip in the 3 2 representation. The invariants in the dimension-5 operator 
are as in case 25. We can get a tribimaximal neutrino mixing matrix by choosing the following 
different options for the vev's: 

i) 






where b), c) and k) will contribute, 
ii) 

/! 

(x) = x, (k) = ui 1 

V o 

leading to the same result as in i). 
iii) 

(n) = Vl ^ j , (k) = Ul ^ 1 j , (C) = u 2 ^ 1 | (103) 

with contributions as in i) plus extra terms from a) and h). 
iv) 

' \ / 1 




:/»; - r 2 | j , («)=ui ^ 1 J . (0 !(2 I J I (1-01) 

with contributions from b), c), i) and k). 
v) 

(p) = v 2 ^ 1 j , («) = ui ^ 1 | , (C) = u 2 | I (105) 

with terms contributing from b), c), i) and k). 

The following corresponding neutrino mass matrices will be possible: 




i , n 



iii) 



iv) 



f2u\ f[ 9 UiU 2 /( 9 tilW2 

M u = \ f[ gUl U2 hu 2 2 f A u 2 2 ) (106) 

f{gUlU 2 f±u\ fou\ 
H u \ + fl4VlU2 f[ 9 UlU2 f[gUlU2 

M v = \ f' l%Ul U2 hu\ hv\ + f 4 u 2 2 \ (107) 

/( 9 ^2 fivj + f A ul fou\ 

hu\ + fl7V2Ul f[ 9 UlU2 A9U1U2 

M v = ( f[ 9 u lU2 f 2 u\ + hu 2 2 f' 4 u 2 2 + f[ 5 v 2Ul 1 (108) 

f[ 9 uiU2 fiuj + f[ 5 v 2 ui f 2 u\ + foul 
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v) 



M ^ = I fi5 v 2Ui f2u{ + fnv 2 ui f[ 9 uiu 2 I (109) 



/i 5 v 2^1 /i 9 u l u 2 /2«? + fvtVlUl 

all with the same form as in case 23: 




M v = \ (3 a 5 \ (110) 



This will again give tribimaximal mixing, if <5 = 7 + /3 — a is fulfilled. 

The selected cases are chosen, because there exists no other result than these. Only the two 
possible neutrino mass matrices extracted here, can be obtained in this group. We can therefore 
cover all the cases of the group, by just studying these two forms of the neutrino mass matrix, 
with the ij) in either a 3i or a 32 representation. 



10 Conclusion 

The Frobenius group, with its four 3 dimensional irreducible representations, seems like a suitable 
group to use, in order to study the family symmetry of neutrinos and charged leptons. 

The approach we have used here is similar to the approach used to study the tetrahedral 
group in [7]. It is interesting to notice how similar the two groups are with similar irreducible 
representations, Kronecker products and Clebsch-Gordon decompositions of these, as well as 
the mass matrices extracted from the Lagrangian. However, the results gained for the vacuum 
alignments are very different indeed, and the match between the neutrino and charged lepton 
sector seems to be improved for this Frobenius group. However, to find out which vev is the 
right to use, we need to place further constraints on the model. 

In this paper, we have not addressed the sequestering problem. When the Higgs of the 
dimension-4 operator acquires a vev in the (1, 1, 1) T direction, the T13 symmetry is broken down 
to Z3. As the possible Higgs fields of the dimension-5 operators acquire vev's, some also break 
T13 down to Z3, but in every single case there is at least one vev that breaks it down to Z 2 - The 
Z 2 subgroup does not commute with the Z3 subgroup, and the vacuum alignments can not be 
sequestered from each other, unless the interaction terms in the Higgs potential, connecting the 
fields that break the group down differently, vanishes |18j . The tribimaximal mixing is therefore 
a lowest order approach, in which we need to study allowed deviations within the T13 structure. 
This will impose extra constraints on the vacuum alignments of the Higgs fields, changing the 
mixing matrix, which thereby might also change the fit between the mass matrices gained from 
the neutrino and charged lepton sector, respectively. 

Certain articles have reached modifications to the tribimaximal mixing [19J [20J. These could 
be used for this theory, to find a more constrained match between the two sectors. Also, we 
have in this paper assumed CP conservation. A violation of this could be considered with an 
appropriate mixing matrix. 

Furthermore, we have assumed a specific vacuum alignment of the Higgs field in the dimension- 
4 operator. Other vacuum alignments could be studied. 
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Note : After our first submission, we have come to learn, that another paper, [22J, has car- 
ried out studies of the group X13 as a model for tribimaximal mixing, using a supersymmetric 
approach. 
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